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1 . INTRODUCTION 


It has been discovered experimentally that a so-called ” plateau" 
region exists in the vortex system which trails from a lifting wing. 

The decay of the vortex due to viscous or turbulent shear is very slow 
in the plateau so that the maximum tangential speed in the vortices 
remains very nearly constant for some distance downstream of roll-up and 
then begins to decrease, becoming inversely proportional to the square 
root of the distance far downstream. 

It is possible that the delay in vortex decay, exemplified by the 

plateau region, is due, at least in part, to nonequilibrium turbulence 

in which the magnitude of turbulent shear stress takes a significant 

period of time to catch up with turbulent energy. It is shown here, 

however, that the existence of the plateau region can be explained merely 

by the necessity for the vortex tangential speed profile to change from 

an inviscid profile (inviscid except for the core) to the similarity 

profile far downstream. The inviscid profile is a function of span 

2 3 

loading calculated a la Betz 1 , so the extent of the plateau region is 
also a function of span loading. 

It is assumed in the following analysis that the flow in the trailing 
vortex is approximately analogous to the time-dependent flow of an infinitely 
long vortex, so the radial and axial momentum equations are not used. This 
assumption usually results in good qualitative agreement in the near field 
just downstream of roll-up and should be close to a three-dimensional 
solution in the far field flow. Solutions have been obtained for both 
constant and variable eddy viscosity models. 
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2. VORTEX DECAY - CONSTANT EDDY VISCOSITY 

The equation for an axisymmetric, infinite vortex assuming a constant 
eddy viscosity is 

£-%($-*!?) <■> 

To facilitate solution of this equation, the Independent variables are 
changed : 


t = t/t 

o 

N = r 2 /4v T t 


and the equation becomes 


lY _ 
5r 


N 


hx _ 

3N 


= N 


dN 2 


( 2 ) 


2 

subject to the boundary conditions y(l,N o ) = Y^(r , y(t,o) sa 0, and 

Y(t,°°) = Y 0 « This is, of course, a linear differential equation 
solvable by separation of variables: 


Y= G(N)T(t) 


dT , . 2_ 

T — + \ T = 0 


( 3 ) 


N + N + \ 2 G = 0 
dN 2 dN 


where \ represents one or more constants 


( 4 ) 
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The solution to Eq. (3) is 


T = 




( 5 ) 


and a solution to Eq. (4), which is a special case of Rummer* s equation , 
is given by 


G = 


\ 2 -l 


dN 


X 2 -l 


/ X 2 -N. 
(N e ) 


( 6 ) 


Thus, a solution of Eq. (2) which satisfies the boundary conditions is 


1 


a> 




,X 2 -1 


dN 




(N 



A A 

-N 1 -N 2 2 -N 

= 1 - e + — Ne + “ (2N - N^)e + etc. 

T i 

T 


( 7 ) 


The first two terms of this equation represent Lamb's solution for decay 
of an infinite line vortex. 

The relationship between vortex radius and circulation for elliptic 

2 

loading is given according to the Betz theory by 

r/b = { [tt/8 - (1/4) sin 1 (1 - y ) l ^ 2 ]/'y] - j (1 - 

( 8 ) 


The variable N at completion of roll-up (t= 1) is an unknown constant 
2 

times (r/b) . Thus, when fitting the solution Eq. (7) to the inviscid 
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Eq. (8), the coefficient of N in the exponent must be solved for as well 
as the coefficients A 9 . Since the resulting set of algebraic equations 

r 

is nonlinear, the procedure for solving for the unknown coefficients is 
necessarily iterative. An example of the solution Eq. (7) for elliptic 
loading is 


= i + -1.6(10 r/b ) 2 / t | . 1.5804(10 r/b ) 2 

Y ' " 2 

o i t 


|~ 7 .4890 ^ 10 r/b ^ 2 _ 5,9912 ^ 10 r/b ^ 4 j 


^ -15.1114 ^ 10 r/b y 


24.1782 


^ 10 r/b ^ 4 _ 6,4475 ^ 10 t/b j j 
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, T 3 \ T / ' T 2 l " ) - V - J 


3.8600 ^ 


+ [^8^2540 (10_ 

r/b \ 2 58.4129 /10 r/b\ 4 

) T 3 I ^ j 

46.7302 

— i 

/l0 r/bV* 

12.4614 { 

/l0 r/b\ 

6 0.9969 /10 r/b\ 10 l 

2 

T 

l T ) 

^ T ( 

, - ) 

\ T / J 

j“ 7.4438 

(“) 

2 29.7751 j 

4 ! 

T 

flO r/b \ 

l - ) 

4 ( 31.7601 ^10 r/b^ 6 

12.7035 

/10 r/b\ 8 

. 2.0326 /l0 r/b\ 10 

0.1084 / 10 r/b\ 12 1 ) 

2 

T 

l T ) 

T V 

r I 

\ T ) J I 


( 9 ) 
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This equation fits the inviscid solution fairly well, as shown in Fig. 1, 

except for the viscous core region. That the actual experimental profile 

corresponds to the inviscid Betz solution (except for the viscous core) 

has been shown to be true by Rossow^. 

A previous solution for a constant viscosity, non-potential vortex 

6 

was obtained by Kirde . He obtained a solution involving confluent 

hypergeometric functions for which the initial tangential velocity varies 
- 1/2 

3s r , as does the inviscid Betz solution for small radius for elliptic- 

loading. However, Kirde 1 s solution is not valid for large radius as the 

1/2 

circulation in his initial vortex varies as r for large r rather than 
approaching a constant value as does the solution Eq. (9). 

The tangential velocity profiles in the vortex for various values of 
time following roll-up (At = t - 1) are shown in Fig. 2. For large values 
of time At, the profile approaches the Lamb solution. The maximum tangential 
speed as a function of time (or downstream distance) is shown in Fig. 3, 
along with experimental data from water channel, wind tunnel, and flight 
tests. For low Reynolds number, the distance data has been modified by 
a function of Reynolds number, as described in Ref. 7. Although the 
constant viscosity assumption results in a profile that does not match 
experiment far downstream, still the maximum velocity curve seems to match 
the data very well, showing both the near field plateau and the similarity 
decay region far downstream. Thus, although nonequilibrium turbulence may 
have a hand in shaping the plateau region, it has been shown here that the 
existence of the plateau can be attributed at least primarily to the 
adjustment of the vortex profile from the inviscid to the viscous similar 
solution. 

Solutions for constant eddy viscosity were also obtained for two other 
span loadings. The circulation profile for triangular loading, for which 



CIRCULATION ~ y/y, 



RADIUS ~r/b 


Fig. 1. Initial circulation profile, elliptic loading 
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Fig. 2. Tangential speed profiles as a function of time, elliptic loading. 
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Fig, 3. Maximum tangential speed vs, downstream distance, elliptic loading 




9 


the inviscid solution is given by 


y" = 4 r/b, r/b <0.25 


( 10 ) 


is shown in Fig. 4. The solution of Eq. (2) in this case is 


y = 1 + e 
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The velocity profiles corresponding to Eq. (11) are shown in Fig. 5, 
a span loading such that 


Y- [1 - (2y/b) 2 ] 3 , 


For 


( 12 ) 


CIRCULATION ~y/ 



Fig. 4. Initial circulation profile, triangular loading 




TANGENTIAL SPEED RATIO ~V/V 



Fig. 5 . Tangential speed profiles as a function of time, triangular loading 
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Y = 1 + e 
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(13) 


The circulation and velocity profiles for this case are depicted in 
Figs. 6 and 7, respectively. 

The effect of span loading on the initial maximum tangential velocity 
level and duration of the plateau is shown in Fig. 8. The triangular 
loading results in a maximum tangential speed immediately after roll-up 
which is less than half* that for elliptic loading. This result has also 
been shown experimentally 1 . Thus, the hazard to trailing aircraft would 
be reduced for generating aircraft which exhibit this type of span loading. 
Far downstream, the solutions merge and the hazard far downstream, which 
is less than in the plateau, becomes independent of span loading. 




Fig# 6* Initial circulation profile, intermediate loading. 



TANGENTIAL SPEED ~ V/V 



Fig. 7. Tangential speed profiles as a function of time, intermediate loading. 
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Fig. 8. Maximum tangential speed vs. downstream distance, loading comparison. 
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3. VORTEX DECAY - VARIABLE EDDY VISCOSITY 

As in Ref. 7, in order to simulate the variation of turbulent 
shear stress with vortex radius, the eddy viscosity is approximated, 
using the mixing length concept, as 

v t - a2 r2 l r f; I 

With this eddy viscosity model and changing variables from t, r to 

2 2 

t - an d N=r/4a y t, Eq. (3) of Ref. 7 becomes 

T S = { 4 I N !sf- Y ' + ‘' /a \) N ~2 + N In < 15 > 

dN 

Equation (15) was solved using a Crank-Nicholson finite difference 
scheme. The velocity profiles for a vortex trailing from an elliptically 
loaded wing for various times following roll-up are shown in Fig. 9. 

Figure 10 illustrates the maximum tangential speed in the vortex as a 
function of time or downstream distance. Because of the fact that the 
initial profile [Eq. (8)] was completely inviscid in this case, the 
plateau does not appear, and the initial decay in the vortex is approxi- 

-1/3 g 

mately AX , substantiating C. E. Brown’s theory for elliptic span 
load profiles. Far downstream, however, the effect of span load gradually 
disappears, and the maximum tangential speed eventually decays as 
as shown. 

Also shown in Fig. 10 are the results of a calculation in which the 
initial inviscid profile was modified to incorporate solid body rotation 



TANGENTIAL SPEED ~ y/VNt 



Fig* 9* Tangential speed profiles as a function of time, elliptic loading, variable eddy 
viscosity* 
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o 


Fig. 10. Maximum tangential speed vs. downstream distance, elliptic loading, variable eddy 
viscosity. 
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in the core (illustrated in Fig. 11). For this example, the core radius 
chosen iis quite small so the extent of the plateau is short. The end of 
the plateau region is illustrated by the dotted line in Fig. 10. Experi- 
mental data for approximately elliptic span loads (Fig. 3) show relatively 

larger core sizes, a longer plateau, and a very much shorter region in 

-1/3 

which the maximum speed varies as 

Velocity profiles for a triangular span loading [Eq. (10)] are 

illustrated for various distances downstream in Fig. 12. The inviscid 

Betz solution for this case results in a tangential speed which is constant 

from the center of the vortex out to one-fourth the span. Viscosity and 

turbulent shear quickly round the corners of the initial profile, but the 

maximum tangential speed stays nearly constant out to At - 0.5, as shown 

- 1/2 

in Fig. 13 by the extensive plateau. Similarity and the £X decay 
follow the plateau very quickly, in contrast to Fig. 10 and in harmony 
with the data of Ref. 1. 

Figures 14 and 15 exhibit the core radius and circulation values as 

functions of downstream distance for the three calculated cases. The 

existence and relative extent of the plateau is also clearly shown by 

the elliptic viscous core and triangular loading curves for core radius 

in Fig. 14. The core radius for the inviscid elliptic case initially 
0 6 

grows as AX * , which is slightly different from Brown’s expectation of 
2/3 

AX . The difference is probably due to the fact that the profile shape 

continually changes and therefore does not satisfy the criterion of 

1/2 

similarity. Far downstream the radius grows as . The difference 

in core radius circulation between elliptic and triangular span loading 
as evidenced in Fig. 15 is striking. The circulation ratio, r^/T , for 
the elliptic case would increase in the plateau region with increase in 




Fig- 11. Initial tangential 



TANGENTIAL SPEED ~ y/Vn^ 



VORTEX RADIUS ~yNr 


Fig. 12. Tangential speed profiles as a function of time, triangular loading, variable eddy 
viscosity. 



MAXIMUM TANGENTIAL SPEED 


10 



Fig* 13# Maximum tangential speed vs, downstream distance , triangular loading, variable 
eddy viscosity. 



100 



Fig. 14. Vortex core radius vs. downstream distance, loading comparison, variable eddy viscosity 




VORTEX CORE CIRCULATION 



Fig. 15. Core radius circulation vs. downstream distance, loading comparison, variable eddy 
viscosity. 
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initial core size, but would usually be less than the far downstream 
value of 0,41 which results from the turbulent eddy viscosity model 
used here. The effect of the span loading in the near field and of 
variable eddy viscosity in the far field is to keep the value of core 
radius circulation far below that predicted by the Lamb solution for 

Q 

a potential vortex (f^/T = 0.715). Therefore, previous reports of 

puzzlingly low values of circulation ratio are easily explained. 

4. CONCLUSIONS 

The existence of a plateau region in which the trailing vortices 
exhibit a very slow decay of maximum tangential speed has been verified 
by the use of both constant and variable eddy viscosity models. While 
the shape of the plateau is probably influenced by the effects of axial 
flow and nonequilibrium turbulence, it has been shown that the existence 
of the plateau can be explained by the presence of the viscous core and 
by the necessity for the initial, nearly inviscid tangential velocity 
profile to undergo transformation to a turbulent similarity profile far 
downstream. 

Variation in span loading is shown to affect the initial maximum 
tangential speed in the vortex as well as the duration of the plateau. 
Triangular span loading would reduce the hazard to trailing aircraft from 
that due to more conventional loading since the initial tangential speed 
is much lower. This can be explained by the fact that the inviscid profile 
for triangular loading exhibits a finite tangential speed at the vortex 
centerline and an extensive plateau results even without an initial 
solid body core. Elliptic span loading, however, results in an inviscid 
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profile having an infinite speed at the centerline of the vortex, and 
a viscous core is thus necessary for the existence of a plateau. 
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APPENDIX 

In order to match the variable eddy viscosity profiles with experiment, 

it will be necessary to adjust coefficients used in the calculation of 

2 

initial profiles. For elliptic loading, the relationship between 
circulation y and independent variable N is 


N = K l{ [ 4 ' 2 si“’ X Vl - Y 2 ] / Y - fVl - Y*} (Al) 


where 


r/b - V” 7 


(A 2) 


vb/r = ( Vk77tt)y /Vn? 

O i 


(A3) 


r x 

O o 
U b b 




(A4) 


the initial core size will also need to be adjusted. For triangular loading 


N = K 2 y 2 (r/b < 0.25) 


(A5) 


r/b = (1/4 VK^) V¥^ 

(A6) 

Vb/r o = (2 Vk^tt) y/VTTt 

(A7) 

Ir x 


a \]uJ> ~ * y^T/ 32 K 2 

(A8) 


the values of and K 2 used in the examples were 755 and 3.4346 respectively. 



